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ABSTRACT 

We analyze the resolved stellar populations of 473 massive star-forming galaxies at 0.7 < z < L5, with multi- 
wavelength broadband imaging from CANDELS and Ha surface brightness profiles at the same kiloparsec resolution 
from 3D-HST. Together, this unique data set sheds light on how the assembled stellar mass is distributed within 
galaxies, and where new stars are being formed. We find the Ha morphologies to resemble more closely those 
observed in the ACS I band than in the WFC3 H band, especially for the larger systems. We next derive a novel 
prescription for Ha dust corrections, which accounts for extra extinction toward Hii regions. The prescription 
leads to consistent star formation rate (SFR) estimates and reproduces the observed relation between the Ha /UV 
luminosity ratio and visual extinction, on both a pixel-by-pixel and a galaxy-integrated level. We find the surface 
density of star formation to correlate with the surface density of assembled stellar mass for spatially resolved regions 
within galaxies, akin to the so-called “main sequence of star formation” established on a galaxy-integrated level. 
Deviations from this relation toward lower equivalent widths are found in the inner regions of galaxies. Clumps 
and spiral features, on the other hand, are associated with enhanced Ha equivalent widths, bluer colors, and higher 
specific SFRs compared to the underlying disk. Their Ha /UV luminosity ratio is lower than that of the underlying 
disk, suggesting that the ACS clump selection preferentially picks up those regions of elevated star formation 
activity that are the least obscured by dust. Our analysis emphasizes that monochromatic studies of galaxy structure 
can be severely limited by mass-to-light ratio variations due to dust and spatially inhomogeneous star formation 
histories. 
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1. INTRODUCTION 

Studies of galaxy evolution from the peak of cosmic star 

formation to the present day have matured tremendously over 

the past two decades. Initially, efforts focused on the detection 

and selection of distant blobs of light by means of color selec- 

tions (Steidel et al. 1996; Franx et al. 2003; Daddi et al. 2004), 

of which generally brighter subsets were confirmed spectro- 

scopically. Gradually, galaxy-integrated photometry over ever 

increasing wavelength baselines, from the rest- frame UV to 

infrared (IR), enabled inferences of the global stellar popula- 

tion properties such as stellar mass, age, and star formation 
rate (SFR; Papovich et al. 2001; Forster Schreiber et al. 2004; 

Shapley et al. 2005). Alongside developments on the stellar pop- 
ulation front, deep-field observations carried out with the Hubble 
Space Telescope ( HST) have played a pivotal role in opening up 
the distant universe for resolved, morphological analyses (e.g., 
Williams et al. 1996; Giavalisco et al. 2004; Rix et al. 2004; 


Koekemoer et al. 2007). Initially, such look-back surveys were 
predominantly monochromatic, or at most spanning a limited 
number of bands probing the rest-frame UV emission, which is 
dominated by young O and B stars. 

The first HST - based extensions into the rest-frame optical 
of high-redshift galaxies already date back more than a decade 
(Thompson et al. 1999; Dickinson et al. 2000), but until recently 
were limited to small sample sizes due to the small field of 
view of the NICMOS camera. The installment of the Wide 
Field Camera 3 (WFC3) during Servicing Mission 4 changed 
the situation dramatically. Benefiting from WFC3’s enhanced 
sensitivity and larger field of view, the Cosmic Assembly 
Near-infrared Deep Extragalactic Legacy Survey (CANDELS; 
Grogin et al. 2011; Koekemoer et al. 2011) has been mapping 
the rest-optical properties of galaxies since the peak of cosmic 
star formation at high resolution in five disjoint fields on 
the sky (GOODS-South, GOODS-North, COSMOS, UDS, 
EGS). Exploiting initial subsets of the CANDELS legacy data 
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set, several works have contrasted monochromatic structural 
measurements to galaxy-integrated stellar population properties 
(Wuyts et al. 2011b; Weinzirl et al. 2011; Szomoru et al. 
2011; Bell et al. 2012; Wang et al. 2012; Kartaltepe et al. 
2012; Bruce et al. 2012). Overall, these authors found strong 
correlations between structure and stellar populations (i.e., a 
“Hubble sequence”) to be present out to at least z ~ 2.5 (see 
also Franx et al. 2008; Kriek et al. 2009a). Star-forming galaxies 
(SFGs) residing on the SFR-mass “main sequence” relation 
(Noeske et al. 2007; Elbaz et al. 2007, 2011; Daddi et al. 2007) 
tend to be the largest at their mass, and are best characterized by 
exponential disk profiles, while quiescent galaxies lying below 
the main sequence have higher Sersic indices, i.e., are more 
bulge-dominated (Wuyts et al. 2011b; although see Barro et al. 
2013 for interesting sub-populations deviating from this overall 
trend). This implies bulge growth and quenching are intimately 
connected. 

The molecular gas content of normal SFGs is known to in- 
crease rapidly with look-back time, to gas mass fractions of 
-0.33 (0.47) at z ~ 1.2 (2.2) (Tacconi et al. 2010, 2013; 
Daddi et al. 2010). Simple stability arguments predict such gas- 
rich systems to be prone to gravitational collapse on scales of 
— 1 kpc. Indeed, such features, and more generally irregular mor- 
phologies, have been abundantly reported in high-redshift SFGs 
(albeit on monochromatic, and often rest-UV or Ha images; see, 
e.g., Elmegreen et al. 2004, 2009b; Genzel et al. 2008, 2011). 
Analytic work as well as hydrodynamic simulations have pro- 
posed that the giant clumps may migrate inward via dynamical 
friction and tidal torques, thereby providing an alternative chan- 
nel to bulge formation to the conventional merger scenario (e.g., 
Noguchi 1999; Immeli et al. 2004a, 2004b; Bournaud et al. 2007; 
Dekel et al. 2009). In part, the efficiency of such violent disk 
instabilities in driving significant bulge growth is subject to the 
longevity of the clumps, which may be limited by vigorous out- 
flows driven by internal star formation (e.g., Newman et al. 2012; 
Genel et al. 2012; Hopkins et al. 2012). However, several studies 
argued that even if the clumps do not remain bound, torques in 
unstable disks will still lead to an enhanced gas inflow rate with 
respect to stable configurations (Krumholz & Burkert 2010; 
Bournaud et al. 2011; Genel et al. 2012; Cacciato et al. 2012). 

Additional clues on the emergence of bulges and the role or 
nature of clumps can come from studies of the resolved stellar 
populations within galaxies. Here too, Advanced Camera for 
Surveys (ACS) + WFC3 imaging from CANDELS and pre- 
existing campaigns offer the capability to expand sensitively on 
pioneering work by, e.g., Abraham et al. (1999), Elmegreen 
et al. (2009a), and Forster Schreiber et al. (2011a, 2011b). 
Wuyts et al. (2012) carried out seven-band stellar population 
modeling on a (binned) pixel-by-pixel basis for 649 massive 
(logM* > 10) SFGs at 0.5 < z < 2.5, the largest and 
only mass-complete sample of SFGs subject to such a detailed 
analysis to date. Translating the internal variations in intensity 
and color to spatial distributions of more physically relevant 
quantities such as stellar mass, SFR, age, and extinction, these 
authors found high-redshift SFGs to be smoother and more 
compact in mass than in light, with color variations driven by 
a combination of radial extinction gradients and spatial (short- 
term) fluctuations in the star formation history (see also Forster 
Schreiber et al. 2011a, 2011b; Guo et al. 2012; Lanyon-Foster 
et al. 2012; Szomoru et al. 2013). In particular, Wuyts et al. 
(2012) and Guo et al. (2012) found regions with enhanced 
surface brightness with respect to the underlying disk to be 
characterized by enhanced levels of star formation and younger 
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ages than interclump regions at the same galactocentric distance. 
Typical inferred clump ages of 100-200 Myr at z ~ 2 imply that 
the clumps correspond to short-lived star-forming phenomena, 
possibly limited in lifetime by stellar feedback. If inward clump 
migration is taking place, this should happen efficiently on 
timescales similar to the orbital timescale. Radial age gradients 
of clumps (Forster Schreiber et al. 2011b; Guo et al. 2012) may 
signal such migrational processes to be at play. At z — 1, which 
is the epoch we focus on in this paper, the physics leading to 
regions of excess surface brightness and locally enhanced star 
formation may be a mix of the above violent disk instabilities 
frequently studied at z — 2 and more conventional processes 
known from nearby SFGs. That is, gas fractions can still be 
sufficiently high to cause gravitational collapse while timescales 
become long enough for other instabilities, such as spiral density 
waves, to arise (the latter are only in rare cases seen at z ~ 2; 
see Law et al. 2012). We note, as we did in our previous work, 
that while for simplicity we occasionally use the same short- 
hand “clump” terminology as Wuyts et al. (2012), the specific 
diagnostic identifying enhanced surface brightness regions does 
not discriminate between round knots or spiral features. 

Despite the richness of the HST broadband data sets, the 
multi-wavelength sampling available at high resolution is still 
modest in comparison to state-of-the-art spectral energy distri- 
butions (SEDs) on a galaxy-integrated level. The latter span up 
to 30+ medium and broad bands from 0.3 to 8 m (e.g., Ilbert 
et al. 2009; Whitaker et al. 2011; Spider et al. 2012), com- 
plemented further by mid-IR-to-far-IR photometry, and probe 
dust-obscured star formation within at least the more massive 
SFGs (e.g., Magnelli et al. 2013). The performance and lim- 
itations of broadband SED modeling as a tool to infer stellar 
population properties have been abundantly tested on synthetic 
observations of simulated galaxies (Wuyts et al. 2009; Lee et al. 
2009; Pforr et al. 2012; Mitchell et al. 2013). It is clear that, 
while four to seven HST broad bands may be sufficient to re- 
construct stellar mass distributions, 13 some level of degeneracy 
(e.g., between age and dust extinction) will be inherent to the 
inferences made on a pixel-by-pixel basis when SEDs span the 
relatively narrow wavelength range from observed B to H band. 
Additional empirical constraints would allow us to test and build 
on the more model- sensitive findings from Wuyts et al. (2012). 

The 3D-HST legacy program (van Dokkum et al. 2011; 
Brammer et al. 2012) and GOODS-North grism survey GO- 
11 600 (PI: B. Weiner), together covering all five CANDELS 
fields, provide two such empirical constraints. First, the grism 
data yields spectroscopic redshifts for thousands of galaxies, 
eliminating an important source of uncertainty that propagates in 
all derived stellar population properties. Second, each dispersed 
galaxy image can be considered as a continuum with superposed 
resolved line emission maps at HST resolution. After subtracting 
the continuum underlying the Ha emission of z ~ 1 galaxies, we 
are therefore left with an independent SFR diagnostic probing 
the same kiloparsec scales as the CANDELS multi-wavelength 
broadband imaging. The additional information furthermore 
allows us to better constrain the effects of extinction. Nelson 
et al. (2012, 2013) demonstrated the power of this technique by 
contrasting the stacked one-dimensional Ha and H - band surface 
brightness profiles. They found the extent of the ionized gas and 


13 Stellar mass estimates to an accuracy of ~0.2 dex can arguably even be 
obtained on the basis of two broad bands, provided they bracket the 
Balmer/4000 A break (Bell & de Jong 2001; Taylor et al. 2011; and see 
Forster Schreiber et al. 2011a for an application of this technique to the 
resolved mass distribution in a z ~ 2 galaxy). 
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stellar light profiles to be similar for small galaxies. The ionized 
gas component of the larger galaxies, on the other hand, exhibit 
larger scale lengths than the stellar component, consistent with 
an inside-out growth scenario. As is the case for rest-optical 
light profiles, the Ha emission of SFGs in their sample is best 
characterized by an exponential disk profile. 

In this paper, we build on the analyses by Nelson et al. 
(2012, 2013) and Wuyts et al. (2012) by combining the re- 
solved Ha information on z ~ 1 SFGs with pixel-by -pixel stel- 
lar population modeling of their multi-wavelength broadband 
photometry. The present work has strong connections to a very 
comprehensive exploration of Ha and SFR profiles throughout 
the two-dimensional (2D) SFR-mass space (E. Nelson et al. 
in preparation). While restricting the analysis to systems with 
well-detected line emission, we expand the sample with respect 
to previous work by exploiting the data in ah five C ANDELS / 
3D-HST fields (see the sample description in Section 2). We 
discuss the methodology to derive stellar population properties 
such as star formation and surface mass density in Section 3. 
Here, we place a special emphasis on dust corrections, and in 
particular extra extinction toward the H ii regions from which 
the Ha line emission originates, using UV+IR based SFRs as an 
additional calibrator. After a visual impression of characteristic 
features in the Ha and broadband images (Section 4), we quan- 
tify the correspondence between Ha and broadband morpholo- 
gies at different wavelengths (Section 5.1), compare emission 
line diagnostics to stellar population properties derived from 
broadband information alone (Section 5.2), and finally contrast 
the Ha and corresponding star formation properties of galaxy 
centers, clumps/ spiral arms, and underlying disks within our 
sample of SFGs with log M * > 10 (Section 5.3). 

Throughout this paper, we quote magnitudes in the AB sys- 
tem, assume a Chabrier (2003) initial mass function (IMF) and 
adopt the following cosmological parameters: (O m , Qa, h) = 
(0.3, 0.7, 0.7). 

2. OBSERVATIONS AND SAMPLE 

HST broadband imaging from CANDELS and pre-existing 
surveys, together with HST grism spectroscopy from 3D-HST, 
form the core data set on which this paper is based. In addition, 
our sample definition makes use of the wealth of ancillary 
data collected in the five CANDELS/3D-HST fields, including 
multiple medium-band and broadband imaging campaigns from 
the ground, and space-based photometry from Spitzer/ IRAC, 
Spitzer/ MIPS, and Herschel/ PACS. We derived the galaxy- 
integrated stellar masses and SFRs based thereupon following 
identical procedures to Wuyts et al. (2011b). That is, the SFRs 
that enter the sample selection criteria are based on a ladder of 
SFR indicators, using UV+PACS for PACS-detected galaxies 
and UV+MIPS 24 /x m for MIPS-detected galaxies; SFRs from 
SED modeling (which also yield the stellar masses) are used 
for sources without IR detection. The 3D-HST catalogs with 
consistent multi-wavelength photometry that serve as input to 
the galaxy-integrated SED modeling are presented by R. Skelton 
et al. (in preparation). 

2.7. CANDELS HST Imaging 

CANDELS provides deep WFC3 imaging of five disjoint 
fields on the sky: GOODS-South, GOODS-North, COSMOS, 
UDS, and EGS, totaling approximately 800 arcmin 2 . Point- 
source limiting depths vary from H = 21.0 mag in CANDELS/ 
Wide to H = 27.7 mag in CANDELS /Deep (the central halves 
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of the GOODS fields). We refer the reader to Grogin et al. 
(2011) for an overview of the survey layout and Koekemoer 
et al. (2011) for details on the data reduction. Other HST 
imaging that enters our analysis includes the GOODS ACS 
campaigns (Giavalisco et al. 2004) and subsequent ACS epochs 
taken as part of the HST supemovae search, COSMOS/ACS 
(Koekemoer et al. 2007), AEGIS/ACS (Davis et al. 2007), and 
ERS/WFC3 (PI: O’Connell) observations. Together, these yield 
£435V606*775Z85o7i25#i60 photometry in the GOODS fields, 
with an additional To98 coverage of the ERS region (top third 
of GOODS-South), and T105 coverage of the CANDELS/ 
Deep regions. For the remaining fields, SEDs are sampled by 
V 606 huJnsHm photometry. 

2.2. 3D -HST Grism Spectroscopy 

Brammer et al. (2012) describe in detail the specifications of 
the 3D-HST legacy program. Briefly, 3D-HST 14 covers three 
quarters (625 arcmin 2 ) of the CANDELS treasury survey area, 
and naturally our sample will be drawn from the region where 
the two data sets overlap. In this paper, we make use of the 
two orbit depth WFC3/G141 grism exposures, but note that 
in parallel two to four orbits of ACS/G800L grism data were 
taken as part of the program. The WFC3 /G141 grism data reach 
a 5a emission-line sensitivity of ~5 x 10“ 17 erg s -1 cm -2 and 
cover a wavelength range from 1.1 to 1.6 /xm. We make use 
of the v2. 1 internal release by the 3D-HST team, which differs 
most significantly from the data handling described by Brammer 
et al. (2012) in that the half-pixel dithered exposures are 
combined by interlacing rather than drizzling, thereby reducing 
the correlated noise. The extraction of grism spectra followed 
the steps outlined by Brammer et al. (2012). As part of the 
extraction pipeline, a spectral model is constructed that, once 
convolved with the F MOW pre-image of the galaxy that serves as 
a template, best fits the observed grism spectrum. In this process, 
neighbors are treated simultaneously to reduce contamination 
effects. Redshifts are fitted to the combination of available 
broadband and grism information. In practice, for the sample 
of line-emitting galaxies analyzed in this paper, the ancillary 
broadband data only serves in the redshift determination to 
prevent line misidentihcations. 

2.3. Sample Definition 

The adopted redshift range of 0.7 < z < 1.5 for our sample 
is dictated by the requirement that the Ha emission falls within 
the WFC3/G141 wavelength coverage. We next apply the same 
basic criteria to select massive SFGs as Wuyts et al. (2012), 
namely 1 ogM* >10 and a specific SFR (SSFR) > 1/ bubble • 
The depth of the H- band selected parent catalogs guarantees 
completeness down to this mass limit (and in fact more than 
an order of magnitude below). Over the entire CANDELS /3D- 
HST area with coverage by the G1 41 grism spectroscopy and at 
least four HST broad bands, this amounts to 1844 massive z ~ 1 
SFGs. 

Subsequent selection criteria serve to optimize the quality 
of our inferences on resolved Ha properties at the expense of 
compromising the mass completeness of our SFG sample. For 
our final sample, we require galaxies to have secure redshifts 
with well-covered Ha emission detected at the 8a significance 
level. We furthermore apply a conservative screening against 


14 The G141 grism coverage of the GOODS-North field from program 
GO- 11 600 (PI: B. Weiner) is incorporated into 3D-HST as the observational 
strategy is nearly identical to that of 3D-HST. 
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Table 1 

Sample Selection of Massive SFGs at 0.7 < z < 1.5 with Ha Maps from 3D-HST 


Selection 

N (Total) 

N (GOODS-S) 

N (GOODS-N) 

N (EGS) 

N (UDS) 

N (COSMOS) 

0.7 < z < 1.5 & log M* > 10 

3115 

609 

656 

340 

692 

818 

... & SSFR >1/ ^Hubble 

1844 

359 

381 

206 

433 

465 

. . . & Ha S/N > 8 

636 

109 

176 

96 

130 

125 

... & no contamination or continuum residuals 

506 

95 

143 

72 

110 

86 

... & Final sample a 

473 

85 

129 

69 

108 

82 


Note. a Compact or heavily point-source-dominated sources with an X-ray counterpart are excluded from the final sample. 


objects whose grism spectra are contaminated by neighboring 
sources, as identified based on products of the 3D-HST pipeline, 
supplemented with a visual inspection. In this step, we also weed 
out a small number of sources with prominent residuals over 
the full wavelength range after subtraction of the continuum 
model produced by the 3D-HST extraction and fitting pipeline. 
Finally, we exclude 33 bright X-ray detected objects for which 
the compact nature of their HST imaging could hint at a 
significant non- stellar contribution to both the broadband and 
the Ha line emission. Before the latter cut, the overall fraction of 
X-ray detected sources in the Ha sample exceeds the respective 
X-ray detected fraction determined for the complete parent 
sample of massive SFGs by a factor of 1.5 (~12% versus 
~8% 15 ). A possible explanation for this difference may be a 
tendency for actively star-forming, line-emitting galaxies to 
more frequently host active galactic nucleus (AGN) activity 
(e.g., Santini et al. 2012; Trump et al. 2013). In addition, since 
the Ha and [N ii] emission lines are blended in the WFC3 grism 
data (see Section 3.2), the enhanced [Nii]/Hq' ratios of AGN- 
hosting galaxies may push more of them over the signal-to-noise 
threshold for line emission. In our final sample, the only X-ray 
detected sources remaining are located in GOODS-South (16) 
and GOODS-North (7). We note that in these deep X-ray fields 
(4 Ms and 2 Ms exposure, respectively), not all X-ray detected 
sources are necessarily AGNs. Moreover, no evidence for point 
source contamination is seen in their HST imaging. We break 
down our final sample and the intermediate selection steps by 
field in Table 1. 

It is important to note that our final sample of 473 massive 
SFGs spans the entire range in mass and SFR of typical main- 
sequence galaxies above log(M*) = 10 (see Figure 1). The inset 
panels in Figure 1 present a closer look at the relation between 
the Ha sample and the underlying complete parent population 
of massive SFGs at 0.7 < z < 1.5. At first glance, the samples 
are well matched in redshift, mass, star formation activity, 
obscuration (SFRir/SFRuvX axial ratio, and size, in terms of 
dynamic range spanned as well as the distribution within that 
range. In more detail, the Ha sample shows subtle biases against 
the lowest SSFR systems, against the most obscured galaxies 
at a given SFR, and, related, against the most edge-on galaxies. 
The similarity between the two samples is encouraging, as it 
implies that the analysis of the Ha sample presented in this 
paper can reveal generic insights for the entire massive SFG 
population during the transition from the peak of cosmic star 
formation to the initial phases of its decline. 


15 The percentages quoted are for the sum of all five CANDELS/3D-HST 
fields. While the absolute percentages are field-dependent due to varying depth 
of the X-ray imaging, the relative factor of 1.5 between the Ha and parent 
sample is similar for all fields. 



Figure 1. Location of our massive SFG sample (red) in SFR-mass space, 
overplotted on the distribution of the overall galaxy population at 0.7 < z < 1.5. 
Crosses mark X-ray-detected sources. Inset panels show the distribution in 
redshift, stellar mass ([M 0 ]), specific SFR ([yr — 1 ]), and obscuration level 
(SFRir/SFRuv) for the Ha sample analyzed in this paper (red) compared 
to the complete sample of SFGs above log M* = 10 (black). 

(A color version of this figure is available in the online journal.) 

3. METHODOLOGY 
3.1. Modeling of the Broadband SEDs 

The three basic steps toward resolved SED modeling are 
point-spread function (PSF) matching, pixel binning, and stellar 
population synthesis modeling of the individual spatial bins. 
Each of these steps is described in depth by Wuyts et al. (2012). 
Briefly, we work at the WFC3 Higo resolution of O'/ 18 and used 
the Iraf PSFMATCH algorithm to build kernels to bring all 
shorter wavelength images to the same PSF width and shape. 
We next applied the Voronoi binning scheme by Cappellari & 
Copin (2003) to group neighboring pixels together in bins so as 
to achieve a minimum signal-to-noise level of 10 per bin. We 
applied the Voronoi binning scheme to the WFC3 Hi^o maps. 
Photometry in the other (PSF-matched) ACS and WFC3 bands 
was measured within the identical bins of grouped pixels in 
order to guarantee that the colors entering our analysis probe 
consistently the same physical regions. The resulting multi- 
wavelength photometry per spatial bin is then fed to EAZY 
(Brammer et al. 2008) to derive rest- frame photometry and 
FAST (Kriek et al. 2009b) to fit stellar population synthesis 
models from Bruzual & Chariot (2003) with identical settings 
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to Wuyts et al. (2011b, 2012). Since all galaxies in our sample 
have spectroscopically confirmed redshifts from the Ha line 
detection in the grism data (and in half of the cases confirmed 
independently by ground-based spectroscopic campaigns), we 
fix the redshift to its spectroscopically determined value in our 
SED fitting, hence reducing the number of free parameters 
by one. In addition, we assume the stars have a fixed, solar 
metallicity. While this assumption is frequently adopted in the 
literature on stellar populations of massive z ~ 1 galaxies, we 
caution that additional constraints and further investigations are 
needed to address the validity and impact of this assumption 
for resolved studies on kiloparsec scales. We follow Salim 
et al. (2007) in adopting parameter values marginalized over the 
multi-dimensional grid explored in SED fitting, as this approach 
proved most robust (compared to, e.g., adopting the least squares 
solution) in the limit of poorly sampled SEDs. Additional notes 
on the reliability of the broadband SED modeling and its impact 
on the results presented in this paper are discussed in the 
Appendix. 

3.2. Extracting Ha Emission Line Maps 

In order to extract Ha emission line maps, we follow Nelson 
et al. (2012, 2013), and subtract the continuum model from 
the observed 2D grism spectrum. The resulting residual image 
is then mapped to the CANDELS frame using redshift and 
astrometric information, and contains the emission line surface 
brightness distribution without imposed prior on its morphology. 
At the spectral resolution of 8v ~ 1000 km s -1 , Ha and [Nil] 
AA 6548+6583 are blended. As we lack additional constraints, 
we apply throughout the paper a simple downscaling of the 
observed emission line flux by a factor of 1.2, and refer to 
this quantity as the Ha flux. In reality, [Nn]/Ha ratios may 
vary between galaxies as well as spatially within galaxies (see, 
e.g., Liu et al. 2008; Yuan et al. 2012, 2013; Queyrel et al. 
2012; Swinbank et al. 2012; Jones et al. 2013; N. M. Forster 
Schreiber et al. in preparation). However, restricting the above 
spectroscopic samples to the same redshift range as adopted 
in this paper, the scatter in [Nil] /Ha is substantial compared 
to any systematic trend, if present, with galaxy mass above 
log(M*) = 10. Furthermore, the [Nn]/Ha gradients reported 
based on adaptive optics assisted observations are typically 
shallow (N. M. Forster Schreiber et al. in preparation). We 
conclude that a higher order correction than the uniform scaling 
factor we apply is not justified by the present data. 

Again following Nelson et al. (2013), we apply a wedge- 
shaped mask to regions that could potentially be affected by 
[Sn] AA67 16+6731 line emission (redward from Ha along the 
dispersion axis), which can mimic the appearance of an off- 
center clump in the Ha maps. 

3.3. Dust Corrections to the Ha Emission 
3.3.1. Birth Clouds and Diffuse Interstellar Dust 

Proper extinction corrections are critical for the physical 
interpretation of dust- sensitive SFR tracers. This applies to 
short wavelength (rest-UV) broadband indicators, motivating to 
a large extent programs such as CANDELS. While at the longer, 
rest-optical wavelengths, attenuation laws predict a significantly 
suppressed impact by dust; this may not be the case for the Ha 
emission at 6563 A, depending on the geometry of dust and 
stars. Ha emission originates from Hu regions immediately 
surrounding young star- forming regions, which are known to 
be often associated with enhanced levels of obscuring material. 
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As such, the nebular emission emerging near massive O stars 
that have not yet dispersed or escaped from their dust-rich 
birth clouds will be subject to extra extinction with respect to 
continuum light at the same wavelength that is produced by the 
bulk of the stars (no longer embedded in the molecular clouds 
in which they once formed). The observed (i.e., attenuated) Ha 
flux then relates to the intrinsic flux Fn a ^ nt as 

Fua , obs = T Ha , in, KT 04 *- l(r 0 - 4A ™ (1) 

where A cont represents the attenuation by diffuse dust in the 
galaxy and A extra represents the attenuation happening locally 
in the birth cloud. The latter have a negligible covering factor 
and therefore affect no other galaxy light than that emerging 
from the respective star-forming region itself. The geometrical 
picture sketched here relates intimately to the framework of the 
Chariot & Fall (2000) model and subsequent refinements by, 
e.g., Wild et al. (2011), Pacifici et al. (2012), and Chevallard 
et al. (2013). 

Empirical evidence for the need of an extra extinction 
correction to Ha (i.e., A extra ^ 0) was first presented for a sample 
of nearby starburst galaxies by Calzetti et al. (1994, 2000). Also 
at larger look-back times, evidence for differential extinction 
between nebular regions and the bulk of the stars has been 
mounting, from arguments based on the observed Ha equivalent 
widths (EWs; van Dokkum et al. 2004; Forster Schreiber et al. 
2009), comparisons of multi- wavelength SFR indicators (Wuyts 
et al. 2011a; Mancini et al. 2011), and measurements of the 
B aimer decrement (Ly et al. 2012; Price et al. 2013). While 
showing a median consistency with the local calibration by 
Calzetti et al. (2000; A extra = 1.27A con t), the uncertainties and 
sample sizes used in those studies did not allow the authors to 
discriminate between a constant A extra (i.e., constant birth cloud 
properties independent of the optical depth from the diffuse 
component) and one that scales proportionally with the diffuse 
column as A extra a A cont . Conceptually, one can think of the 
diffuse attenuation as being determined by the galaxy’s gas 
fraction, dust-to-gas ratio, large-scale geometry, and orientation. 
The attenuation from birth clouds is expected to share some of 
these dependencies (e.g., dust-to-gas ratio), but not others (e.g., 
orientation, as demonstrated by Wild et al. 2011). If the gas 
fraction of a galaxy increases, all other properties remaining the 
same, nothing will change to A ex tra if this results merely in an 
increased number of star-forming clouds, each with identical 
conditions. If the mass of individual clouds does increase with 
galaxy-integrated gas fraction, it depends on the cloud scaling 
relations whether this impacts A ex t ra (e.g., for a Larson 1981 
scaling relation p c ioud ~ r _1 the optical depth would remain 
unaffected). 

3.3.2. Calibrating Ha Dust Corrections Using 
Multi-wavelength SFR Indicators 

From these considerations, we conclude that the precise func- 
tional form of Aextra may not be a constant or proportional to 
A CO nt? but rather a hybrid form between those two extremes (i.e., 
Agxtra saturating with increasing A cont ). By lack of other direct 
constraints at z ~ 1, we take an empirical approach to cap- 
ture this behavior and define a prescription that yields a small 
scatter and negligible systematic offset between two sets of 
galaxy-integrated SFR measurements. As reference, we adopt 
the ladder of SFR indicators from Wuyts et al. (2011b), which 
is based on UV+IR emission, or SED modeling in cases where 
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Figure 2. Comparison of Ha -based SFR estimates to the reference ladder 
of SFR indicators from Wuyts et al. (2011b), in order of priority based on 
UY + Herschel/ PACS, or Spitzer / MIPS 24 /rm, or U- to-8 /xm SED modeling. 
We mark UV+IR-based SFRs in red and indicate galaxies without IR detections 
with gray circles. Crowded regions of the diagram are displayed with a darker 
hue. Ha SFRs need to be corrected for dust extinction to avoid underestimates 
(top panel). Applying an extinction correction corresponding to what is inferred 
from SED modeling is insufficient (middle panel). Accounting for extra 
extinction toward Hu regions, we find a good correspondence to the reference 
SFRs with modest scatter and without systematic offsets (bottom panel). 

(A color version of this figure is available in the online journal.) 

no IR detection is available (see Section 2). Figure 2 contrasts 
Ha -based SFRs from 3D-HST to the reference SFRs. The top 
panel iterates the obvious necessity for dust corrections to the 
Ha emission. The large sample statistics exploited here also 
demonstrate unambiguously the presence of additional extinc- 
tion toward the nebular regions (middle panel of Figure 2, where 
A e xtra = 0 is assumed and significant systematic underestimates 
for highly star-forming systems are evident). Finally, the bottom 
panel of Figure 2 shows the improved agreement between the 
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SFR measurements once differential extinction is accounted for. 
The adopted prescription, 

A extt a = 0.9A cont -0.15A2 ont , (2) 

yields a scatter of 0.227 dex and a negligible systematic offset 
of —0.035 dex. We note that applying the Calzetti et al. (2000) 
prescription for extra extinction (corresponding to A Hcy = 
Acont/0-44, or equivalently Aextra = 1.27A con t) leads to a slightly 
higher systematic offset of 0.066 dex and produces a larger 
scatter of 0.284 dex. If we were to adopt the Calzetti et al. 
(2000) prescription, all objects with high inferred A v values 
from broadband SED modeling would have dust-corrected Ha 
SFRs systematically in excess of the reference indicator. We 
note that the need for an extra extinction correction is a result 
that is largely driven by the more actively SFGs. For most 
of them, UV+IR (i.e., bolometric, rather than dust corrected) 
measurements of the SFR are available. Those galaxies that lack 
IR detections in the CANDELS/3D-HST fields generally suffer 
less obscuration and are therefore less sensitive to the precise 
dust correction applied (although applying a constant A ex tra> 
even to sources with A con t near zero, would lead to overestimated 
Ha SFRs at the low-SFR end). 

3. 3. 3. Ha / UV Ratios 

We now assess the validity of the calibration to correct Ha 
for dust extinction (Equation (2)) by considering an independent 
dust-sensitive diagnostic, namely the Ha'/UV luminosity ratio. 
Here, we define the UV luminosity as the rest- frame luminosity 
L 2800 = vL y (2800A), which we compute with EAZY. In 
principle, the Lh«/Tuv ratio is not uniquely dependent on 
dust extinction, but may also vary among galaxies that differ 
in their relative number of O and B stars due to differences 
in the star formation history or initial stellar mass function 
(see, e.g., Meurer et al. 2009). The reason is that the Ha 
emission is powered by the ionizing radiation from O stars 
(with typical lifetimes of ~7 Myr), while the rest- frame UV 
light receives its major contributions from both O and B stars 
(i.e., stellar lifetimes up to ~300 Myr). The left-hand panel 
of Figure 3 visualizes the evolution in Lh«/Tuv for a set of 
BC03 stellar populations models with varying star formation 
histories. Here, we calculated the intrinsic Hey luminosities from 
the rate of H ionizing photons in the BC03 models. Applying the 
recombination coefficients for case B from Hummer & Storey 
(1987) for an electron temperature of T e = 10 4 K and density 
of n e = 10 4 cm -3 , this gives 

log(L H Jerg s- 1 ]) = log(V Lyc [s- *]) - 11.87 (3) 

where V Lyc is the production rate of Lyman continuum photons 
from the stars. It is noticeable that, except for very young stellar 
populations or abruptly declining star formation histories, the 
balance between the intrinsic (i.e., unattenuated) Ha and UV 
emission quickly converges to a value of around —1.7. We re- 
mind the reader that the frequently used Kennicutt (1998) SFR 
conversions, which assume a 100 Myr old constant star forma- 
tion population, correspond to the same ratio. In the following, 
we adopt log(LHa/Tuv) = — 1-7 as the expected intrinsic ra- 
tio in the absence of dust. We make the plausible assumption 
that any observed variation in L Ha /Tuv will be dominated by 
extinction effects, and that any (minor) contributions from IMF 
variations or short-term fluctuations in the star formation history 
do not depend systematically on Ay. 
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Figure 3. Left: dependence of the intrinsic (i.e., extinction free) Ha/UV (2800 A) luminosity ratio on age since the onset of star formation for a range of star formation 
histories (exponentially declining or increasing and constant star formation). The Kennicutt (1998) SFR conversions from Ha and UV assume a population with 
constant star formation ongoing for lOOMyr. Right: given a stellar population as assumed by Kennicutt (1998), the observed Ha/UV luminosity ratio depends on the 
visual extinction Ay as indicated by the black lines, which correspond to three different models for extra extinction toward Hu regions: no extra extinction (solid), 
double Calzetti (Ana = A cont /0.44; dashed), and limited extra extinction (Equation (2); dash-dotted). The latter model best describes the locus of observed galaxies 
(small red dots). 


(A color version of this figure is available in the online journal.) 


It is then possible to represent different prescriptions for extra 
extinction toward the nebular regions by lines in a diagram 
of Lft a /Luv versus A v (right-hand panel of Figure 3). In the 
absence of extra extinction, the Hey /UV ratio increases linearly 
with increasing A v , while it remains relatively constant (in 
fact drops slightly) with increasing A v for a so-called “double 
Calzetti” prescription. Neither prescription provides a satisfying 
description of the locus of observed 3D-HST galaxies in the 
diagram (red symbols). Instead, the curve corresponding to 
Equation (2) matches the data well. In Section 5.2, we will 
demonstrate that our prescription for Ha dust corrections, the 
associated dependence of Hey/UV ratios on extinction, and 
in general the above geometrical picture of a two-component 
dust distribution still holds on the kiloparsec scale resolution of 
the HST data. A direct implication of the observed increase in 
Hey/UV ratios with Ay is the fact that Ha luminosities are less 
sensitive to dust than UV luminosities. 

Finally, as an additional verification of the applicability 
of our Hey dust prescription (Equation (2)), we repeated the 
comparison of SFR estimates and our analysis of the L Ha /^uv 
versus Ay diagram for the upper and lower tertiles of the axial 
ratio distribution of our galaxies separately. We confirm that a 
good correspondence remains, with only marginal systematic 
shifts. While the sense of the minor inclination dependence 
qualitatively agrees with the findings by Wild et al. (2011) for 
nearby galaxies (edge-on systems requiring relatively less extra 
extinction because a larger fraction of the obscuring column 
toward nebular regions is contributed by diffuse dust), we opt not 
to introduce an extra inclination term in our prescription, given 
its small amplitude compared to the scatter among galaxies of a 
given axial ratio. 

4. CASE EXAMPLES 

Before quantitatively analyzing the Ha properties of our 
sample of massive z ~ 1 SFGs, it is insightful to visually 
inspect a set of case examples. Figure 4 shows a subset of 12 
SFGs, illustrating the diversity in size, color, and morphology 
represented in the sample. Despite their diverse appearance, 
some key trends that also characterize the sample as a whole 
can be noted. 


First, the central regions of most (and in particular the larger) 
galaxies are more pronounced in the observed H than in the I 
band. This reflects the negative color gradients (red cores) found 
in previous work and more compact stellar mass distributions 
based thereupon (e.g., Szomoru et al. 2011, 2013; Wuyts et al. 
2012; Guo et al. 2012). Second, the I - band images appear 
clumpy or feature spiral signatures more frequently than the 
rest-optical H- band images, echoing findings that the clumps/ 
spiral arms are bluer than the underlying disk (Wuyts et al. 2012; 
Guo et al. 2012). 

New with respect to these previous studies is that we can 
now also contrast the Ha morphologies to the short- and long- 
wavelength broadband morphologies. A visual inspection tells 
us that overall the Ha morphologies resemble more closely the 
I- band images than the H- band surface brightness distributions. 
Both the rest-UV light and Hey emission are dominated by the 
youngest stellar populations, unlike the rest-optical emission 
to which the bulk of the stars (including older populations) 
contribute. It is, however, not a trivial observation because of 
the different wavelengths probed and dust effects discussed 
in Section 3.3. Moreover, exceptions exist in which the Hey 
morphology does not trace the I - band light. For example, the I- 
band emission of GN_9 129 is entirely dominated by a clump that 
is not the most prominent Hey clump in the object. In addition, 
for the smaller galaxies we are hard-pressed to make a statement 
on whether the visual morphological match of the Hey map is 
closer to I than to the H band. Galaxies such as GS_22173 
just appear compact in all bands considered. These results are 
reminiscent of the conclusions drawn by Nelson et al. (2012) on 
the basis of a sample of 56 EW- selected galaxies in 3D-HST. 

5. RESULTS 

5.7. Ha Emission Traces the Rest-UV 
Better Than the Rest-optical Light 

In order to quantify the above visual impressions, we analyze 
the morphological match between Ha and H and Hey and I 
band surface brightness distributions respectively by cross- 
correlating the images. In doing so, we allow minor shifts of 
the galaxy centroid before mapping the grism image to the 
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Figure 4. Gallery of case examples from the massive z ~ 1 SFG sample. PSF-matched three-color postage stamps sized 374 x 374 are composed of 1775/125/660 
for galaxies from the GOODS fields, and 614 625 #160 for galaxies in EGS/UDS/COSMOS. Below, we show the surface brightness distributions in I, H, and Ha, 
respectively (at natural resolution, with a slight smoothing applied to the Ha maps for visualization purposes). Blue, star-forming regions present in the / band generally 
dominate the Ho: emission as well. Central peaks in surface brightness (i.e., “bulges”) appear more prominently in the //band. 

(A color version of this figure is available in the online journal.) 


astrometric frame of the CANDELS broadband images. This 
accounts for uncertainties in the position of the grism spectra 
on the detector and degeneracies between on the one hand 
the systemic redshift and on the other hand morphological k- 
corrections between the Ha? line map and the galaxy’s F140W 
light distribution that was used as template in the redshift 
determination. For each band, we adopted the shift leading 
to the largest correlation coefficient, in most cases limited to 
0-3 pixels. 16 The ratio of I-to-H band correlation coefficients 
then serves as a quantitative measure of how much better the 


16 One pixel corresponds to 23 A in wavelength, and 07 06 spatially. 


Ha morphology corresponds to the I - band light distribution 
compared to the //band. Figure 5 demonstrates what we inferred 
by eye: a majority (65%) of galaxies in our sample shows a 
stronger cross-correlation between the / and Ha morphologies 
than between the H and Ha morphologies. This difference 
is most pronounced among the more extended systems. For 
galaxies with semimajor axis lengths larger than 3kpc (as 
measured with GAFFIT by van der Wei et al. 2012), more 
than 75% shows a better match between Ha and the / band than 
between Ha and the //band. GN_96, GN_373, and GN_4952 
in Figure 4 serve as examples of 3-5 kpc sized sources with 
log(cross J x Ha /cross _H x Ha?) ~ 0.05-0.07. For the largest 
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Figure 5. Top: histogram of the ratio of cross-correlation strengths between 
the Ha and 7- band surface brightness distribution and the Ha- and 77- 
band surface brightness distribution, respectively. The vertical dashed line 
marks the median. The distribution is asymmetric with a tail toward positive 
values of log(cross J/cross_H), representing objects whose Ha profiles greater 
resemble those observed in the 7 band than in the 77 band. Bottom: the closer 
correspondence between the Ha- and 7-band images is most prominent for 
the larger galaxies in our sample (with R e being the effective radius measured 
by fitting a Sersic profile to the 77-band light profile). More compact galaxies 
generally look similar at all considered wavelengths (7, 77, and Ha). 

SFGs at z ~ 1, with R e > 5kpc, this fraction increases 
further to above 90% (see, e.g., EGS_7021 and COS_16954 
in Figure 4). We stress that the above statements on the 
resemblance between line emission maps and broadband images 
of different wavelengths are of a comparative nature. Substantial 
pixel-to-pixel variations exist in the flux ratio between Ha? and 
any single broad band (even at short wavelengths), preventing a 
simple scaling from a monochromatic broadband image to Ha? 
without additional information. In the next section, we explore 
how multi-wavelength broadband information improves this 
situation and furthermore allows us to dust correct Ha maps 
to star formation maps. 

5.2. Comparing Resolved Broadband and 
Emission-line Diagnostics 

5.2.1. Resolved Dust and SFR Calibrations 

We now proceed to translate the direct observables (i.e., sur- 
face brightness profiles in Ha and multiple broad bands) to 
more physically relevant quantities following the resolved SED 
modeling and Ha extinction correction methodologies outlined 
in Section 3. Figure 6 contrasts various Ha? -based diagnostics 
to stellar population properties inferred using broadband infor- 
mation alone. We place each spatial bin of each galaxy in the 
respective diagrams and consider the locus of observed points, 
marked by contours. Due to noise, some individual spatial bins 
have negative Ha fluxes and therefore do not appear in the 
plot. These pixel bins are, however, included in the binned me- 
dian and central 50th percentile statistics marked in red. In a 
resolved analogy to Figure 3, we confirm the distribution of 
observed Ho?/UV ratios to relate to the visual extinction Ay 
inferred from resolved SED modeling as would be expected 
from Equation (2). In contrast, A extra = 0 or A extra = 1.27A cont 
do not reproduce the observed locus of galaxy spatial bins in 
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the diagram. Consequently, it comes as no surprise that, when 
applying dust corrections according to Equation (2) to the Ha 
emission, we find the Ho? -based Xsfr t0 correspond well to 
dust-corrected UV-based estimates. To compute the latter, we 
converted the rest- frame 2800 A luminosity L 28 oo to an unob- 
scured UV SFR based upon the Kennicutt (1998) conversion 
SFRuv, uncorrected = 3.6 x 10 _10 L 2 8oo/^o and finally applied an 
extinction correction based on the broadband SED modeling, 
which assumed a Calzetti et al. (2000) attenuation law. Since 
the visual extinction appears in the quantities on both axes of the 
Ha?- versus UV-based dust-corrected X S fr, we also show that a 
strong correlation remains when plotting the directly observed 
Ha? luminosity surface density (i.e., uncorrected for dust) as a 
function of it would be inferred to be on the basis of the available 
broadband information. 

5.2.2. The Resolved Main Sequence of Star Formation 

Using the analysis of SFR estimates and dust treatment, we 
now turn to Figure 7, which contrasts the Ha? -based SFR surface 
density of individual sub-galactic regions to the surface density 
of stellar mass present at the respective location in the galaxy. As 
such, the diagram can be regarded as a spatially resolved version 
of the SFR-mass plane that has been the focus of many galaxy 
evolution studies in recent years. We find a relation between the 
star formation activity and amount of assembled stellar mass 
to be present for sub-regions within galaxies, as is the case 
for galaxies as integrated units. Similar to the well-established, 
galaxy-integrated main sequence of star formation, we observe 
a near-linear slope: 

log(2 SF R[Mo yr" 1 kpc -2 ]) = -8.4 + 0.95 log(X*[M 0 kpc -2 ]), 

(4) 

as fit to the spatial bins with log(X*) < 8.8. Akin to galaxy- 
integrated observations at the high-mass end by, e.g., Whitaker 
et al. (2012), the “resolved main sequence” shows a tendency to 
flatten at the high stellar surface mass density end. Consequently, 
a linear fit to all spatial bins spanning the entire range in 
X*, yields a slightly shallower slope, of 0.84. We conclude 
that, even if global galaxy properties such as halo mass or 
cosmological accretion rate drive the stellar build-up within 
galaxies, this happens in such a way that on average the ongoing 
and past-integrated star formation activity track each other on 
sub-galactic as well as galactic scales. This is at least the case 
down to the ~ 1 kpc scales probed at the WFC3 resolution, albeit 
with a significant scatter and possible exceptions (i.e., reduced 
specific SFRs) at the highest stellar surface mass density end. 

In the remainder of the paper, we will investigate the deviation 
from a perfect Xsfr cx X* relation more closely and illustrate 
where variations in local specific SFRs occur spatially within 
galaxies. Analogous to the galaxy-integrated main sequence 
of star formation, the clumps/spiral arms featuring excess 
star formation with respect to the smooth underlying disk 
(Section 5.3.1) can be thought of as sub-galactic equivalents 
to the starbursting outliers above the main sequence. 

5.3. Resolved Ha Properties in the 
Two-dimensional Profile Space 

5.3.1. The Nature of Excess Surface Brightness Regions 

With the resolved Ha? properties in hand, we now revisit 
the conclusion by Wuyts et al. (2012), based on broadband 
SED modeling, that clumps in massive high-redshift SFGs 
correspond to star-forming phenomena: regions that exhibit an 
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Figure 6. Comparison of resolved galaxy properties based on Ha vs. resolved stellar population properties derived from broadband information alone. Contours 
indicate the locus of all spatial bins of the galaxies in our sample. Red bars with black dots mark the central 50th percentile and median, respectively. With open blue 
circles, we mark the median results for the subset of galaxies in the GOODS fields, where the high-resolution multi-wavelength sampling is best, with seven bands. 
Black histograms illustrate the distribution of spatial bins. The dust-corrected Ha SFRs of individual spatial bins within galaxies compare favorably to estimates from 
the UV or SED modeling (top panels). Bottom left: a good correlation remains when plotting the directly observed Ha luminosity surface density (i.e., without dust 
correction) as a function of what would be inferred on the basis of broadband information, demonstrating the Esfr agreement is not just due to the appearance of Ay 
in the quantities on both axes. Bottom right: the Ho:/UV luminosity ratio of individual spatial bins within galaxies shows a similar dependence on visual extinction as 
observed in a galaxy-integrated sense in Figure 3. 

(A color version of this figure is available in the online journal.) 


enhanced level of star formation with respect to the underlying 
disk. Given short timescales for these spatial fluctuations in 
local star formation, the same substructure is not necessarily 
present in the distribution of stellar mass as would be expected 
if the observed irregular morphologies were all due to merging 
activity. 

Following Wuyts et al. (2012), we employ co-added nor- 
malized profiles to infer the characteristic star formation prop- 
erties of pixels belonging to the galaxy’s central regions, off- 
center clumps (or spiral arms), and the (interclump) outer disk 
(Figure 8). In each of these diagrams, we consider all spatial 
bins of all 473 SFGs as being part of one meta-galaxy, repre- 
senting the massive SFG sample as an ensemble whose resolved 
stellar populations we then dissect. We apply identical criteria 
to Wuyts et al. (2012) in order to separate spatial bins belonging 
to the three regimes: 

[center] lo g(R/R e ) < —0.5 (5) 

[disk] lo g(R/R e ) > -0.5 & (6) 

logff/S,) < 0.06 - 1.6 log (R/R e ) - log (R/R e ) 2 
[clump] log (R/R e ) > -0.5 & (7) 

log (Z/Z e ) > 0.06 - 1.61og(/?/7? e ) - log (R/R e )\ 


where X/X e is the rest-frame 2800 A surface brightness normal- 
ized by the average surface brightness level within the effective 
radius R e . As in Wuyts et al. (2012), we define radii on ellip- 
tical apertures whose center and axial ratio are derived from 
the stellar mass maps. Since the distance from the center of a 
spatial bin to the center of mass of the galaxy can in principle 
be arbitrarily small, we note that the precise galactocentric radii 
at small 1 og(R/R e ) —0.5 becomes meaningless. Throughout 
the paper, we therefore treat all pixels with lo g(R/R e ) < —0.5 
as belonging to the inner regions of galaxies without further 
regard to their precise value of R/R e . In the top panels of 
Figure 8, we first demonstrate that the conclusions drawn by 
Wuyts et al. (2012) based on broadband information alone also 
apply to our Ha sample at 0.7 < z < 1.5. Off-center clumps 
identified at short wavelengths by an excess in surface brightness 
are characterized by bluer rest-optical colors than the underlying 
disk. Taking one step further away from the direct observables, 
stellar population synthesis modeling of the multi-wavelength 
HST colors implies that they correspond to sites of enhanced 
star formation activity (SSFR) with respect to the underlying 
disk. 

The Ha EW can be considered as a complementary, less 
model-dependent diagnostic of the star formation activity. It 
is defined as the ratio of the star formation sensitive line 
emission strength over the continuum flux, which probes the 
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Figure 7. Equivalent to the so-called “main sequence of star formation” for 
galaxies; resolved subregions within galaxies exhibit a correlation between the 
level of ongoing star formation and the surface density of assembled stellar mass 
(plot style as in Figure 6). For reference, the dashed line marks the slope of the 
galaxy-integrated main sequence from Whitaker et al. (2012). 

(A color version of this figure is available in the online journal.) 


bulk of the stars present in the galaxy. As the emission line 
and continuum are measured at the same wavelength and both 
components are subject to the attenuating column of diffuse 
dust in the interstellar medium of the galaxy, extinction effects 
due to the latter divide out. We observe the Ha EW to be 
higher in the clumps/spiral arm regions (with a median rest- 
frame EWhc* ~ 100 A) than in the underlying smooth disk 
(EWh« ~ 70 A), while the lowest EWs are measured at small 
radii. The latter are associated with the regions of highest stellar 
surface mass density and are related to the flattening slope of 
the resolved main sequence at high X* (see Section 5.2.2). The 
implications of the radial trends will be discussed in depth in a 
forthcoming comprehensive overview of Ha and star formation 
profiles as a function of galaxy mass and SFR by E. Nelson 
et al. (in preparation). Here, we consider in more detail the 
characteristic trends with surface brightness at a fixed radius. As 
discussed in Section 3.3.2, an additional extinction correction 
to account for the obscuring material in birth clouds from which 
the Ha line emission originates is required for the EW to be 
an unbiased proxy of the star formation activity. In the bottom 
right panel of Figure 8, we apply this correction and find an 
excess in the clumps’ intrinsic Ha EW by 25% in the median 
with respect to the underlying disk (i.e., interclump regions) at 
the same galactocentric radii. We note that the median excess 
between the EW of clump and disk pixels does not reflect the 
full dynamic range in rest-2800 A surface brightness variations 
at a given radius and the extent of the associated trend in the Ha 
EW. The reason is that the distribution of disk and especially 
clump pixels is dominated in numbers by pixels that lie close to 
the dividing line between the two regimes (see contours). The 
color coding in Figure 8 illustrates that excess EWs by factors of 
two to three are commonly observed when considering the full 
range of surface brightness levels observed at a given radius. 

In order to address whether the clump and spiral arm regions 
selected through our methodology differ in other aspects than 


star formation activity from the diffuse outer disk, we now inves- 
tigate the relation between the Ha and rest-2800 A luminosity 
surface densities (both uncorrected for extinction; Figure 9, top 
panel). As already hinted at by the similar morphological ap- 
pearance of our galaxies in Ha and ACS imaging (Sections 4 
and 5.1), the UV and Ha SFR tracers correlate strongly, both 
spatially and in amplitude. Their measured luminosity surface 
densities prior to any dust correction, however, do not follow 
a relation with a unity slope. Instead, the formal best-fit linear 
relation is 

log(X H a [Lq kpc -2 ]) = 0.2 + 0.8 log(L 2800 [L e kpc -2 ]). (8) 

Two potential contributors to the deviation from a unity slope, 
scatter around the relation, and, in more general terms, behind 
variations in Ha /UV luminosity ratios, are dust extinction and 
variations in the star formation history (see Section 3.3). While 
both Ha and UV light are tracing young stellar populations, 
they probe recent star formation over different timescales 
(around lOMyr and 100 Myr, respectively). The UV and Ha 
radiation are furthermore emitted at different wavelengths, 
and, moreover, continuum and nebular emission may travel 
different paths through the galaxies’ dusty interstellar medium. 
The combination of these effects leads to systematically larger 
Ha/UV ratios with increasing extinction compared to the 
intrinsic (i.e., dust- free) ratio. 

Evaluating the spatial variation of Ha /UV ratios within our 
z ~ 1 SFGs (Figure 9, bottom panel), we find the clumps/spiral 
features to correspond to lower Lh^/Tuv than the underlying 
disk at the same radius. Three effects can potentially contribute 
to this trend. The first relates to measurement uncertainties 
rather than physics. If the true Ha/UV ratio were constant 
across galaxies, measurement errors would preferentially push 
those spatial bins over the clump selection threshold for which 
X 2800 is overestimated. Since the 2800 A luminosity appears 
in the denominator, the measured Ha /UV ratio for those bins 
would consequently be lower than the true value. While possibly 
contributing in a minor way, this effect is unlikely to be by itself 
responsible for the observed trend, given the small uncertainty 
on X 2 8oo ( see Appendix) and the fact that the clumps are also 
confirmed to be physically distinct from the underlying disk 
in quantities that are measured independently of the 2800 A 
luminosity (e.g., the Ha EW). 

A more physical interpretation is that the excess surface 
brightness regions selected from rest-UV ACS imaging are 
not only regions characterized by enhanced star formation, but 
preferentially those local enhancements in star formation that 
are the least obscured. Speculatively, outflows driven by the 
enhanced star formation density in “clumps” (see, e.g., Newman 
et al. 2012) may be responsible for the clearing of gas and dust 
toward these sightlines. 17 

Finally, we cannot exclude that recent variations in the star 
formation history contribute to the lower Ha/UV ratios of 
clumps as well. We note, however, that a recent upturn of 
the SFR in spatial bins that are classified as clumps or spiral 
features would lead to elevated Ha/UV ratios (i.e., relatively 
more O stars compared to the number of B stars), the opposite 
of what is observed. Interpreting the Ha /UV ratios in terms of 
changes in the star formation history would therefore invoke 
a rapidly declining SFR (see Figure 3) and consequently a 


17 See also Calzetti et al. (2005) for a discussion on the increased efficiency of 
dust and gas clearing from high star formation density regions, such as 
starbursts, compared to normal disk galaxies. 
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Figure 8. Co-added normalized rest-UV surface brightness profile of the SFGs in our sample. The vertical dashed line indicates the typical resolution. Contours 
denote the density of spatial bins. The top panels are color-coded by purely broadband-based quantities, namely the rest-frame optical color (U — V) res t (top left) 
and the specific SFR inferred from resolved stellar population modeling (top right). The color-coding in the bottom panels marks the Ha rest-frame equivalent width 
as observed directly (bottom left) and after correction for extra extinction toward the nebular regions (bottom right). Solid lines separate the central, outer disk, and 
clump/spiral arm regimes. Spatial bins in the clump/spiral arm regime are characterized by bluer broadband colors and higher Ha EWs, consistent with the elevated 
SSFR with respect to the underlying disk inferred from the resolved SED modeling. 

(A color version of this figure is available in the online journal.) 


prior (possibly very short-lived) phase of even more intense 
star formation. The latter phase necessarily would have to be 
embedded, since otherwise it would itself be picked up by 
the UV clump selection. As such, both physical explanations 
could go hand in hand, in a scenario where outflows quickly 
clear the gas and dust from clumps, simultaneously reducing 
their star formation activity and obscuration, and consequently 
lowering the observed Ha /UV ratio. If true, this would imply 
that resolved observations at longer wavelengths with ALMA 
or Plateau de Bure Interferometer (PdBI)/NOEMA may reveal 
new/different clumps that correspond to an earlier, embedded 
evolutionary stage of star- forming regions. 


5.3.2. Where Do the Stars Form in z ~ 1 SFGs? 

While the analysis in Section 5.3.1 sheds light on how star 
formation activity (SSFR or EW) varies depending on the 
location within galaxies, it does not reveal how much of the SFR 
integrated over all galaxies is contributed by the center, clump, 
and disk regimes, respectively. Figure 10 serves to address this 
question. We first consider the integrated Ha luminosity of all 
galaxies in the sample, irrespective of whether they have spatial 
bins in the clump regime or not. With 8.3%/79.1%/12.6% of 
the integrated Ha luminosity budget being contributed by the 
central/disk/clump regions (as defined by Equations (5)-(7)), 
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Figure 9. While the Ha luminosity surface density correlates with the rest-frame 
2800 A luminosity surface density of individual bins (top panel), variations in 
the Ha /UV luminosity ratio are present and show a systematic dependence on 
the spatial location of pixels within the galaxies (bottom panel; identical in style 
to Figure 8). The reduced Ha/UV luminosity ratios of clumps with respect 
to the fainter interclump regions may arise from spatial variations in effective 
attenuation (i.e., reduced dust extinction toward those sites of enhanced star 
formation where the rest-UV surface brightness is observed to peak). 

(A color version of this figure is available in the online journal.) 


it is clear that the vast majority of the Ha emission arises 
from the disk component. 18 A similar conclusion can be drawn 
from the breakdown of the rest- frame 2800 A luminosity budget 
with modest differences that can be understood from the trends 
described in Section 5.3.1 and specifically Figure 9. The clump 
pixels account for 18% of the integrated L 2800 luminosity. It 
should be noted that the precise percentages quoted above 
depend significantly on the adopted dividing lines between the 
regions (Equations (5)-(7)), and the fractional contribution of 


18 Here, we accounted for the masked pixels that could potentially be affected 
by [S 11 ] line emission (see Section 3.2), to avoid a biased weighting of 
different radial bins. 
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excess surface brightness regions should be regarded as coming 
from bright clumps and spiral features that show sufficient 
contrast with respect to the extended smooth disk to be detected 
at the kiloparsec resolution of WFC3. For example, if we were 
to lower the dividing line between the disk and clump regions 
by 0.2 dex, bringing it close to the peak of the contours in 
Figure 10, the fractional contribution of clumps in this adjusted 
definition would be 30% to the Ha luminosity and 23% to the 
Ha derived SFR. The latter percentages can be interpreted as a 
rough estimate of the upper limit on the fractional contribution 
by UV bright clumps. 

As discussed extensively in Section 3.3.2, the Ha luminosity 
is not a direct tracer of the total on-going star formation since 
part of the ionized gas emission is filtered out by the intervening 
obscuring material. As discussed by Wuyts et al. (2012), SFGs 
feature radial gradients in dust extinction. This results in an 
enhanced contribution of the inner (log(R/R e ) < —0.5) regions 
of the galaxies to the integrated SFR budget of our sample: 
12.9% when computing the budget based on dust corrected Ha 
emission. A similar contribution of 14.8% is obtained from 
broadband SED modeling. Wuyts et al. (2012) found similar 
fractional contributions for their mass-complete sample of z ~ 1 
SFGs, confirming that our Ha sample is not subject to significant 
biases in this regard. Using broadband HST imaging, Wuyts et al. 
(2012) also extended the analysis to z = 2.5, finding that the 
fractional contribution of clumps to the integrated SFR of SFGs 
increases to ~20% at z ~ 2. 

We conclude that most stars being formed during and since the 
peak of cosmic star formation are born in a disk component. This 
result is consistent with the analysis of SFR profiles presented 
for a smaller 3D-HST sample by Nelson et al. (2013) and is in 
line with a larger body of work based on HST imaging (e.g., 
Forster Schreiber et al. 201 la, 201 lb; Wuyts et al. 201 lb, 2012; 
Guo et al. 2012) and ground-based kinematic measurements 
(e.g., Genzel et al. 2006; Forster Schreiber et al. 2009; N. M. 
Forster Schreiber et al. in preparation; Lemoine-Busserolle & 
Lamareille 2010; Wisnioski et al. 2011; Epinat et al. 2012). 

6. SUMMARY 

We combined CANDELS high-resolution multi-wavelength 
imaging on kiloparsec scales of a sample of 473 SFGs at 
0.7 < z < 1 .5 with maps of the Ha line emission from 3D-HST. 
Together, the CANDELS + 3D-HST data shed light on where 
within galaxies new stars are being formed and in which amounts 
stellar mass has been assembled at each location throughout their 
history prior to observation. As such, the present study builds 
on earlier work from 3D-HST (Nelson et al. 2012, 2013) and 
CANDELS (Wuyts et al. 2012). Compared to the former study, 
we expand the sample to a more representative subset of all 
massive (log(M*) > 10) SFGs at z ~ 1, while compared to 
the latter, we extract sources from all five CANDELS/3D-HST 
fields instead of one. 

In order to optimally exploit the additional constraints from 
Ha on the resolved star formation activity, we started by 
calibrating empirically the dust corrections to the ionized gas 
emission required to match the integrated reference SFRs 
(Wuyts et al. 2011b; for most of our sample UV+PACS or 
UV+MIPS 24 /I m, with values from U- to- 8 /i m SED modeling 
for cases without IR detection). Our empirical calibration 
(Equation (2)) highlights the need for an extra extinction 
correction to Ha compared to the underlying continuum light 
emitted at the same wavelength, although the correction is 
less severe than proposed by Calzetti et al. (2000) based on 
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Figure 10. Co-added normalized rest-UV surface brightness profile of the SFGs in our sample (identical in style to Figure 8). The color-coding of each bin in the 
diagram indicates the percentage of the total integrated Ha luminosity (left panel) or dust-corrected SFR based thereupon (right panel) summed up from all massive 
SFGs in our sample, originating from the respective region in 2D profile space (i.e., galactocentric radius and rest-2800 A surface brightness level). Differences 
between the fractional contribution to the cosmically integrated Ha luminosity and SFR arise from inhomogeneous dust distributions within galaxies: thicker dust 
columns toward the center and slightly reduced obscuration of the rest-UV-selected clumps. The vast majority of newborn stars in z ~ 1 galaxies are being formed in 
the disk. 


(A color version of this figure is available in the online journal.) 


nearby starburst galaxies. Our findings fit in the context of a 
two-component geometrical dust model in which the ionized gas 
emission surrounding young star-forming regions is obscured 
by dust from the birth cloud itself in addition to the column of 
diffuse material spread throughout the galaxy that also affects 
the bulk of galaxy light emitted by older stars (see also Chariot 
& Fall 2000; Wild et al. 2011; Pacifici et al. 2012; Chevallard 
et al. 2013). As an independent test, we confirmed that this 
prescription for extra extinction toward H ii regions reproduces 
the relation between the observed Ha /UV luminosity ratio and 
the effective visual extinction inferred from broadband SED 
modeling. The fact that the observed Ha /UV luminosity ratio 
rises with increasing visual extinction directly implies that Ha 
is a less dust- sensitive SFR tracer than UV emission. 

From our resolved stellar population analysis, we draw the 
following conclusions. 

1. The Ha morphologies of z ~ 1 SFGs resemble more 
closely the ACS I - band morphologies than those observed 
in the WFC3 H band, particularly in the case of large 
galaxies. 

2. The Ha and rest-frame 2800 A luminosities correlate on a 
pixel-by-pixel basis, and after proper dust corrections are 
applied, the SFR estimates based thereupon are internally 
consistent. The Ho'/UV luminosity ratios of individual 
spatial bins also relate to the visual extinction inferred from 
multi-band HST photometry in a manner that is consistent 
with the applied correction for extra extinction toward the 
Hu regions. 

3. We find evidence for the existence of a resolved main se- 
quence of star formation: the rate of ongoing star formation 
per unit area tracks the amount of stellar mass assembled 
over the same area. Its near-linear slope is consistent with 
the one measured for the galaxy-integrated main sequence, 
including an apparent flattening at the high-mass density 
end, associated with the low-EW inner regions of the most 
massive galaxies. 


4. Off-center clumps (or spiral features that show up similarly 
as regions with excess surface brightness) are character- 
ized by enhanced Ha EWs, bluer broadband colors, and 
correspondingly higher SSFRs than the underlying disk, 
implying that they are a star formation phenomenon. Phys- 
ically, they may correspond to regions with elevated gas 
fractions and/or star formation efficiencies (Tacconi et al. 
2013), where part of the global gas disk collapsed through 
gravitational instabilities (see, e.g., Bournaud et al. 2008; 
Genzel et al. 2008; Dekel et al. 2009), although at z ~ 1 sec- 
ular processes on longer timescales, such as spiral density 
waves, are also likely to contribute. In an integrated sense, 
however, the contribution of excess surface brightness re- 
gions to the total amount of star formation taking place 
within z ~ 1 SFGs is limited to 10%-15%, depending on 
the tracer. Most stars being formed between z = 1.5 and 
z = 0.7 are formed in the disk component (see also Forster 
Schreiber et al. 2009; N. M. Forster Schreiber et al. in prepa- 
ration; Nelson et al. 2013). Wuyts et al. (2012) demonstrated 
that this conclusion holds out to z ~ 2, where the fraction 
of total star formation in off-center clumps increases, but 
does not exceed ~20%. 

5 . Aside from spatial fluctuations in the star formation activity, 
an inhomogeneous dust distribution can lead to spatial 
mass-to-light ratio variations, and consequently a range in 
surface brightness levels at a given radius, even for smooth 
stellar mass distributions. The visual extinction gradient 
inferred from broadband SED modeling is to first order 
radial. We furthermore find clumps /spiral features selected 
from ACS imaging to exhibit lower Ha/UV luminosity 
ratios than the underlying disk. Interpreting the Ha/UV 
ratio as a tracer of extinction, this implies that the clumps/ 
spiral features selected in our study are not just regions of 
increased star formation activity, but preferentially those 
observed through a smaller column of obscuring material. 

In the past decade, deep HST imaging has transformed our 
view of high-redshift galaxies from a zoo of faint blobs of light 


14 


The Astrophysical Journal, 779:135 (16pp), 2013 December 20 

into resolved species whose structure and stellar population 
content exhibit strong correlations throughout cosmic time 
(e.g., Franx et al. 2008; Williams et al. 2009; Wuyts et al. 
2011b). More recently, multi-band HST imaging has enabled 
us to translate the distribution of light within galaxies to more 
physically relevant distributions of, e.g., stellar mass or star 
formation (e.g., Wuyts et al. 2012; Guo et al. 2012; Szomoru 
et al. 2013). In this paper, we demonstrated that a more complete 
census of the distribution of baryons within distant galaxies, 
including the ionized gas as traced by Ha (see also Nelson 
et al. 2012, 2013), helps to constrain spatial variations in star 
formation activity and dust extinction. Despite the significant 
strides forward, especially through the synergy between the 
CANDELS and 3D-HST legacy programs, a significant fraction 
of the total baryonic mass budget remains missing from these 
resolved studies. 

In a galaxy-integrated sense, the PHIBSS survey carried out 
with the PdBI revealed that cold molecular gas traced by CO 
accounts for a third of the baryonic mass within z ~ 1.2 SFGs, 
with molecular gas mass fractions increasing further out to 
higher redshifts (Tacconi et al. 2010, 2013; see also Daddi 
et al. 2010). Knowing how the molecular gas is distributed 
within galaxies is crucial to interpret the continuing growth 
of structure: by quantifying the star formation efficiency in the 
central, outer disk, and clump/spiral arm regimes separately, 
identifying the regions with enhanced columns of obscuring 
material, and finally contrasting the total (gas + stellar) baryonic 
mass profiles to the available constraints on enclosed mass from 
kinematics. 

Exploiting the first resolved data set of a distant galaxy 
consisting of three-dimensional (3D) CO and Ha data cubes 
as well as multi-band HST imaging, Genzel et al. (2013) 
demonstrates the power of such a synergy. These authors 
find a correlation between the column of CO line emission 
and the visual extinction inferred from the HST colors, a 
near-linear resolved Kennicutt-Schmidt relation, and consistent 
kinematics of the ionized and molecular gas components. 
The superior sensitivity provided by ALMA and the NOEMA 
extension to PdBI promises to advance this field even further 
through systematically mapping the spatial variations in CO 
line emission as well as the dust-sensitive mm continuum for 
galaxies spanning a wide dynamic range along and across the 
main sequence of star formation. 
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provided by NASA through a grant from the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Incorporated, under 
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APPENDIX A 

ROBUSTNESS OF THE RESOLVED SED MODELING 

As described in Section 2. 1 , the available high-resolution HST 
imagery in the CANDELS/3D-HST fields varies from seven- 
band BVizYJH photometry in the GOODS fields, contributing 
46% of our sample, to a more modest sampling of the SED, in 
VIJH , for EGS, UDS, and COSMOS (objects with less coverage 
were excluded from the sample). Our choice to also include the 
latter three fields in the analysis is therefore necessarily a trade- 
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off between having the richest possible set of information per 
object, and maximizing the sample statistics. 

We verify the validity of our approach in two ways. First, 
we repeat our entire analysis five times, each time using the 
data set of an individual field separately. We revisit the relation 
between Ha and broadband derived properties (Section 5.2), the 
variation of stellar population properties across the co-added 
rest-UV surface brightness profiles (Section 5.3.1), as well as 
the fractional contributions to the Ha luminosity and integrated 
SFR by the central/disk/clump regions (Section 5.3.2). Doing 
so, we confirm that the same generic behavior is seen in each of 
the fields separately, albeit naturally based on smaller number 
statistics. This test boosts our confidence that the conclusions 
drawn in this paper are not driven by just a limited subset of the 
objects entering our analysis. For example, the relative fractional 
contributions of galaxy centers, disks, and clumps or spiral arms, 
vary by at most a few percent between fields. 

Second, we tested the impact of having four rather than seven 
flux points per spatial bin, by repeating the analysis for the 
GOODS-South field using only the ViJH resolved photometry. 
We found that the best convergence between solutions of four- 
and seven-band fitting is obtained when assigning a weight 
Wi = e~ x i /2 to each grid cell and marginalizing over the multi- 
dimensional parameter space explored in the fitting process (see 
Salim et al. 2007 for an extensive discussion). We therefore 
chose to follow this approach throughout the paper, but note that 
consistent conclusions are reached when using the minimum x 2 
value. Comparing the four- to seven-band fitting, we reproduce 
the overall lessons learned on resolved star formation patterns at 
z ~ 1 presented in this paper: the existence of a resolved main 
sequence of star formation with similar slope as the galaxy- 
integrated version, and excursions from this relation in the form 
of off-center regions with excess star formation, or low EW 
regions in galaxy centers. Quantitatively, the stellar surface 
mass density, rest-frame 2800 A surface brightness, and rest- 
frame (U — V)rest color are least affected by the transition 
from seven to four bands, inducing less than 0.07 dex scatter 
and no systematic offsets. The absence of systematics is most 
crucial, as throughout the paper we never treat individual spatial 
bins separately, but instead always consider statistics averaged 
over large numbers of spatial bins taken from all objects in the 
sample. 

Unsurprisingly, the inferred age, SFR, and visual extinction 
are less well constrained by 4-band SEDs, with scatter in the 
four- to seven-band comparison of 0.18 dex, 0.20 dex, and 0.32 
mag, respectively, and strong (anti-)correlations between the 
offsets. However, also for these parameters we measure only 
minor systematic offsets with amplitudes that are a factor six 
smaller than the observed scatter. This propagates to variations 
in the co-added surface brightness profiles and resolved main 
sequence of comparable or even smaller amplitude than the 
variations seen from field to field. 

We conclude that, to the level of detail presented in this paper, 
the above tests support the robustness of our results. However, 
in order to push the study of resolved stellar populations at high 
redshift further, toward extensive analyses on a per object basis, 
additional means to better disentangle age and dust effects would 
be desirable. Two possible avenues are the inclusion of WFC3 
UVIS data (Teplitz et al. 2013), providing access to resolved 
maps of the UV slope /3, and the upcoming James Webb Space 
Telescope /NWSSdim mission, which will enable the application 
of the UVJ diagnostic technique (Wuyts et al. 2007) on sub- 
galactic scales. 
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